Helicase on DNA: A Phase coexistence based mechanism 
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We propose a phase coexistence based mechanism for activity of helicases, ubiquitous enzymes 
that unwind double stranded DNA. The helicase-DNA complex constitutes a fixed-stretch ensemble 
that entails a coexistence of domains of zipped and unzipped phases of DNA, separated by a domain 
wall. The motor action of the helicase leads to a change in the position of the fixed constraint thereby 
shifting the domain wall on dsDNA. We associate this off-equilibrium domain wall motion with the 
unzipping activity of helicase. We show that this proposal gives a clear and consistent explanation 
of the main observed features of helicases. 



Nucleic acid helicases are defined as enzymes that 
translocate directionally through double stranded nucleic 
acid substrates to catalyze the separation of the comple- 
mentary strands. They facilitate various biological pro- 
cesses such as DNA replication, recombination and re- 
pair, RNA transcription, editing and splicingQ]. There 
are several structural varieties of helicases like monomeric 
(e.g. Per A), dimeric (e.g. Rep), trimeric (e.g. RecBCD), 
tetrameric (e.g. RNA polymerase) or closed hexameric 
(e.g. DnaB), but all use the hydrolysis^ of ATP to ADP 
as the preferred source of energy [lj, |2|, y, Lil LA LI III l| ■ 

Bulk behaviour in solutions like average unwinding 
rates, step size, average number of base-pairs opened per 
helicase, etc, are known for a few helicases like Hepatitis- 
C virus helicase^, PcrA 0, DnaB|l E3 and others. 
Much attention has recently been devoted^, 0, Q 
toward a quantitative characterization of RecBCD en- 
zymes, a multifunctional trimeric protein complex (the 
products of the recB, recC and recD genes 0) that par- 
ticipates in the repair of chromosomal DNA through ho- 
mologous recombination. In bacteria, like Escherichia 
coli, RecBCD is involved, e.g., in protection against 
damages by UV or gamma irradiation, and infection 
by bacteriophages. In all cases, the full functionality 
of RecBCD relies on the helicase and the nuclease ac- 
tions of its subunits. The use of single-molecule, mi- 
cromanipulation tools allowed for monitoring in detail 
translocation^, unwinding 0] an d processivity 

(rate of dissociation) 0] of individual RecBCD enzyme 
molecules on dsDNA. Such experiments elucidated sev- 
eral new aspects of helicase behaviour and showed that 
many properties could be more related to general prin- 
ciples than on specific chemical details. In particular 



In general, some nucleoside triphosphate (NTP) is required, e.g. 
T7 gp4 can use GTP, SV40 large T antigen uses non-ATP nu- 
cleotides for unwinding RNA. 9] 



it has been observed that (a) RecBCD unwinds dsDNA 
at a uniform rate, over a wide range of ATP concentra- 
tions, as it moves on one strand, (b) the nuclease activity 
does not affect unzipping, and c) the helicase can work 
in presence of DNA gaps upto certain lengths. More re- 
cently, winding-rewinding for E.Coli Rep helicase-DNA 
complex 8] has been observed at a single molecule level. 

Despite these varieties of experimental findings no 
clear mechanism coupling the motor action and the he- 
licase activity is known yet. To fill this gap, in this pa- 
per we present a simple, but powerful argument, based 
on the principle of phase coexistence [l5l lla. ITtI lisl IT^ | . 
that provides clear and robust explanations to the gross 
observed features. There are a few biological operational 
models built on how a helicase presumably might workQ. 
The analysis reported here gives a thermodynamic basis 
to a model called the "wedge model" according to which 
the motion of the helicase "provides enough force to en- 
able the helicase to destabilize the base pairs at the junc- 
tion by a process resembling the action of a wedge" Q . 
In our proposed mechanism, energy (from ATP) is re- 
quired for translocation activity or the motor action of 
the helicase and not directly for base pair breaking and 
therefore, according to the classification scheme of Rcf. 
@, this corresponds to passive helicases. Additional fea- 
tures required for active helicases are ignored in this first 
study. Our proposal is supported by computer simula- 
tion of an exactly tractable model. To our knowledge, 
this is the first theoretical study of the dynamics of a 
DNA-helicase complex. 

We study the joint dynamics of the helicase and the 
DNA in a two dimensional fork model (Y-Model) [l6lll7| . 
The two strands of DNA are represented by two directed 
and mutually avoiding walks formed by N bases each. In 
two dimensions, on the square lattice (see Fig. ^| this 
implies that the two walks follow the positive direction 
of the diagonal axis (z); in other words the coordinate 
along such direction always increases. The perpendicu- 
lar direction x measures, in unit of the elementary square 
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FIG. 1: (a) A typical configuration of the simulated DNA 
with the helicase (thick rod) as modeled on a square lattice 
(thin lines). Bases are represented by dots and paired bases 
are shown by dotted lines. The position of the fork coincides 
with the last paired base. The arrow is indicating the motion 
of the helicase. Except for the rigid hard-core constraint, 
there is not other direct interaction between the DNA and 
the helicase. (b) Schematic phase diagram in the separation 
(stretch) vs temperature plane. At fixed temperature, a finite 
end-point separation leads to a coexistence of the two phases 
indicated by the filled circles. The interface of the two phases 
is the domain wall. 



diagonal, the distance between the two complementary 
monomers belonging to the two strands. When this dis- 
tance is equal to 1 they are considered in contact: a 
binding energy is gained which is uniform (e = 1) for 
a homogeneous model of DNA (homo-DNA) but chosen 
randomly from two different values (€1,62) for a hetero- 
geneous DNA (hetero-DNA). Notice that due to the geo- 
metrical properties of the lattice the two complementary 
monomers are labeled by the same z-coordinate, as one 
would require for base pairing in DNA. In the Y-model 
unzipping can occur only processively, e.g. bubbles are 
suppressed along the chain: the only conformations con- 
sidered have the first N — m monomers bounded, whereas 
the remaining m are separated in a Y-like conformation. 
The fact that the Y-model does not allow rejoining of the 
unzipped portion of the dsDNA is similar to the geome- 
try observed in the experiment of Ref ^| (see also Fig. 
4 of Ref. H3]). Also, bubbles are suppressed for DNA at 
temperatures much below its melting temperature T m , 
temperature below which the two strands are zipped. In 
the case of homogeneous interaction the exact phase di- 
agram and other static and dynamical quantities can be 
exactly determined (also in the presence of a stretching 
force) Hi 113 • 

The coarse-grained nature of the model[2l| needs to be 
stressed here. Monomers are to be thought of as groups 
of bases, and ignoring helicity or restricting to two dimen- 
sions are more for simplification of the calculation than 
artifacts® . Such coarse-grained models, or even simpler 
ones, are used in various DNA related problems [Til l24| 



[§] The qualita tive features of the unzipping phase transi tion of Refs. 
1 1 -Tit III! Il7j are observed in more complex models in |22| 
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FIG. 2: Zipping probabilities at two different times (after 
thermal averaging) for the homo-DNA (e = I) and hetero- 
DNA (ei = 1, e 2 = 0.5). The solid lines are fits to Eq. The 
inset shows the equilibrium zipping probablities (no domain 
wall) for the case of e = when the ends at z = 400 are kept 
separated. 

and even in analysis of thermal melting of DNA 23] . The 
spirit behind our approach is that the key element that 
can influence universal behaviour of helicase transloca- 
tion is the competition between a Y-fork conformation 
which can be unzipped by paying energy and the move- 
ment of an opening machine. 

Several studies of theoretical DNA models 0, 0, 0, 
have established the existence of a sharp unzip- 
ping phase transition of a dsDNA at a critical stretching 
force applied at one end on the two strands. This im- 
plies that in the conjugate ensemble of a fixed-separation 
constraint for the two strands, there is a coexistence of 
domains of zipped and unzipped phases. For the fixed- 
force ensemble for a homo-DNA in the large length limit 
(N — > oo), the critical force g c for the zipping- unzipping 
transition, and the end separation under a force g > g c 
are given bv|l6l Il7| 

9c (T) = | cosh-V^-lJ.and M=tanh^, (1) 

where temperature T is measured in units of fce/e, kg 
being the Boltzmann constant. From these equations the 
phase coexistence curve (x)/N-vs-T can be determined 
exactly and it is schematically shown in Fig. Under a 
fixed-distance constraint, represented by the vertical line 
in Fig^>, the DNA chain splits into domains of zipped 
and unzipped phases. The length of the unzipped strand 
can be read off from the upper line of the coexistence 
curve. This fact can be checked independently and we 
have verified it not only for the Y-model but also for 
models that allow bubbles. 

We now make our hypothesis. If we think that the he- 
licase, by virtue of its size larger than the the separation 
of the two DNA strands and the excluded volume inter- 
action, acts as a geometrical separator, the helicase-DNA 
complex constitutes exactly a fixed-stretch ensemble. The 
unwinding activity can then be simply associated with 
the motion of the domain wall which necessarily forms 
and follows the motor action of the helicase. In other 



3 



words, the helicase plays a double role: firstly, its pres- 
ence thermodynamically implies the existence of a do- 
main wall, and secondly, its translocation induces a mo- 
tion of the domain wall to reach its equilibrium position. 
Such an approach puts the primary role on the transloca- 
tion motion. The thermodynamic force, that drives the 
domain wall toward its equilibrium position, provides the 
mechanism for base pair opening. 

To verify our hypothesis we have numerically studied 
the DNA model described above mimicking the effects of 
the helicase with a rod of length I (see Fig. QJi). Dynam- 
ics is introduced by a Monte-Carlo procedure. For the 
DNA, one among the 27V monomers is randomly chosen 
and an attempt is made to modify its position with re- 
spect to all the others (which remain fixed). The move 
can in principle increase or decrease by one unit the dis- 
tance between the strands. The move is accepted ac- 
cording to standard Metropolis rules. A Monte-Carlo 
unit time is defined as 2N single- monomer attempted 
moves. The scaling properties of this DNA dynamics, 
also in the presence of a stretching force, have already 
been determined |16|. The helicase moves forward (mo- 
tor action) by unit step, along the — z direction, on the 
DNA, if it is not hindered by the chain configuration 
(excluded volume interaction). The motion is kept uni- 
directional as found for RecBCD in Ref. E2 ( see below 
also). The motion of the helicase is attempted at every 
Monte Carlo step. In our simulation we start with a con- 
figuration where the helicase is attached to one end of 
a dsDNA ( as in experiments of Ref [H El 03 ) ■ The 
temperature is maintained constant for the DNA to be in 
the zipped phase (below T m ). The overall dynamics is off- 
equilibrium. Notice that there is no specific interaction 
but both DNA and helicase dynamics are constrained by 
the excluded volume. 

Throughout the simulation we monitored the position 
of the last zipped base pair (fork position see Fig. [TJi), 
the zipping probability p(s,t) that a base pair at site 
s at time t is zipped, and the average position of the 
helicase at time t. In most simulations the length I of the 
helicase was taken equal to 6 (though we studied up to 
I = 14) whereas 1000 thermal averages were necessary. 
The length N of the strands were varied according to 
the different experiments ranging up to N = 1000. For 
an analysis of the domain wall, the zipping probability 
p(s,t) can be fitted by a function: 

p(s,t) = (1/2) [l-tanh{(«-«o(t))M*)}] (2) 

where so(t) and w(t) are the position and the width, re- 
spectively, of the domain wall. 

In Fig [21 we plot for two different times the zipping 
probability p(s, t) for the homogeneous and for an het- 
erogeneous case. In all cases the domain wall behaviour 
predicted is well fitted by a tanh profile (Eq. [2J. We 
point out that such domain walls are not found if the 
DNA strands are noninteracting, i. e. e = (see inset of 
Fig. HJ. In Fig. [3|we show the instantaneous positions of 
the helicase and the domain wall (computed through Eq. 



as a function of time. Movement starts at t = when 
the helicase is loaded. The two quantities proceed uni- 
formly and cooperatively through the DNA unwinding it. 
The domain wall evolves toward the equilibrium position 
whenever the helicase motion is stopped and this position 
turns out to be not very far away from the instantaneous 
position. This indicates an adiabatic adjustment of the 
domain wall to the instantaneous position of the mobile 
helicase. We stress that had there been no interaction be- 
tween the DNA and the helicase, the latter, as a phantom 
motor, would have moved with the assigned speed with 
the DNA remaining bound (since we are below T m ). It 
is important to compare this motion with the unzipping 
dynamics in a fixed force ensemble which shows a charac- 
teristic scaling behavior[I(| , namely a nonlinear evolution 
m(t) ~ i 1 / 3 , where m(t) is the number of unzipped bases 
at time t. In contrast, and this is the central point of our 
work, we find that the combined dynamics involving the 
excluded volume interaction between the helicase and the 
DNA (but no external force) leads to the uniform motion 
of both the helicase and the domain wall. Our results of 
Fig. should be compared with Fig. 2 of Ref [l3j]. The 
effective velocity is smaller than the unhindered one and 
depends on the size I of the helicase, the temperature and 
the sequence. The I and T dependence can be estimated 
if we think that to move the helicase we have to unzip a 
base pair at the fork position, and to allow the formed 
kink to reach the position of the helicase. The base pair 
is broken with a probability exp(— y) and the kink needs 
a time of order I to reach the helicase. Therefore the ve- 
locity is proportional to [exp(— (confirmed by the 
fit of our data, see the inset of Fig|3]). A strong depen- 
dence on T has also been found in Ref. EH ■ In our sim- 
ulation no sequence dependent nonuniformity in motion 
was ever discernible so long as the heterogeneity was un- 
corrected. The velocity we observed, originating solely 
from the motor action, is a lower bound because any pe- 
riodic conformational change of the helicase E2 during 
its motion (ignored here mainly to illustrate the role of 
the domain wall) would assist the motion of the domain 
wall itself. 

Our proposal lends itself to several predictions. We 
discuss a few here. Since the position so(t) of the wall is 
determined by the location of the stretching constraint 
put up by the helicase, the unzipped part beyond the 
helicase should not affect the action of the helicase. We 
verified this explicitly in our simulation: a part of the un- 
zipped DNA beyond the helicase is chopped off (which 
mimics RecD activity E2) & t arbitrarily chosen times and 
there is neither any change in the nature of the wall (e.g. 
p(s, t)) nor in the rate of unzipping. This agrees with the 
observations in Refs. 0, Il3j that the nuclease activity 
does not modify the unwinding action of RecBCD. Next, 
the width of the wall, as defined by Eq. [21 gives a length- 
scale for the helicase activity. The sequence randomness 
in a heterogeneous chain does not matter if there is no 
correlation beyond this scale, as found both in simula- 
tion and real experiments. We have seen a periodic mod- 
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FIG. 3: Positions of the helicase and the domain wall as a 
function of time for (a,b) homo- and (c,d) hetero- DNA. The 
size of the helicase is I = 6. When the helicase is stopped 
at time= 1000 (indicated by the arrows), the domain wall 
evolves toward the equilibrium position. The inset shows the 
speed of the helicase for different temperatures and lengths, 
with the solid lines representing v — 1.28[exp(— l/T)]/l. 



ulation (not shown) in the helicase motion if the DNA 
sequence is periodic of two pairing energies with period- 
icity larger than the width of the domain wall. Another 
consequence of this scale is that a nick or break smaller 
than the domain wall width will not be recognized, pro- 
viding an interpretation of the observations in Ref 01 
regarding the size of the gap a helicase can negotiate. 

We have also simulated cases where the helicase un- 
dergoes a biased random walk-type motion as expected 
for RecQ.[25j. We introduced a probability P(t) that the 
helicase could step along the — z direction on the DNA 
(when this is sterically acceptable) but also a probabil- 
ity 1 — Pit) that it could move in the opposite direction 
ii. e., away from the domain wall). When P(t) is kept 
fixed at a value Po in the interval 1 > Po > | (Po = 1 
being the case studied in the first part of the paper) we 
find again a uniform motion with velocity related to Pq. 
At Po = \ (random walk) situation changes: unzipping 
does not proceed and the helicase dissociates from the 



double strand. This behaviour is illustrated in 0] where 
we are plotting the position of the helicase as a func- 
tion of time: Pit) decreases linearly from 1 down to | 
and then remains constant. As expected the unwinding 
proceeds non-uniformly until the random walk regime is 
reached. At that time the helicase is discarded and the 
DNA zips again. It suggests, though a bit speculative, 
that a probability affecting the forward motion of the 
helicase could be a phenomenological characterization of 
the gradual wobbliness of the helicase on the track. 

In conclusion, we have shown that associating the 
helicase activity with the domain wall motion in a 
fixed-stretch ensemble accounts for several observed fea- 
tures, as e.g., the uniformity of unzipping, no sequence- 
dependent nonuniformity, and the insensitivity to the 
nuclease action, without any requirement of extra spe- 
cific bond -cutting chemical processes. A domain wall 
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FIG. 4: Positions helicase (I = 6) for an homogeneous DNA 
with a biased probability Pit). Probability decreases linearly 
from 1 for t = to | for t ~ 2500 after which it remains 
constant at 1/2 



also gives a quantitative meaning to the Y-fork in the 
terminology of DNA replication. The underlying ther- 
modynamic basis gives a robustness to the mechanism 
that it could be at work for hexameric helicases also. 
FS was supported by MIUR-COFIN01. FS thanks kind 
hospitality of Institute of Physics, Bhubaneswar. 
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